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Abstract 

We explore the phase diagram of the five-dimensional anisotropic Abelian Higgs model by Monte Carlo 
simulations. In particular, we study the transition between the confining phase and the four dimensional 
layered Higgs phase. We find that, in a certain region of the lattice parameter space, this transition can 
be first order and that each layer moves into the Higgs phase independently of the others ( decoupling 
of layers). As the Higgs couplings vary, we find, using mean field techniques, that this transition may 
probably become second order. 
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1 Introduction 



Coupling anisotropies in gauge theories may lead to fundamental changes in their phase diagrams. Although 
this yields problematic theories if applied to four-dimensional models, due to breakdown of Lorentz invari- 
ance, higher-dimensional models with anisotropic couplings may give rise to theories of physical interest. 
This programme originated in 1984 by Fu and Nielsen |||] who considered a five-dimensional pure U(l) gauge 
theory on the lattice. The main idea is that, for certain values of the couplings for the n extra dimensions, 
four-dimensional layers may be formed within the (4 -|- n)— dimensional space. The corresponding phase is 
called layered phase and one of its main characteristics is that exhibits confinement in the extra dimensions. 
The U(l) higher-dimensional model has already been studied to some extent with lattice techniques (||2|, 
[^), leading to the establishment of the existence of a layered Coulomb phase. | 

The possibility that physical space-time is not four-dimensional has been broadly referred to in the 
bibliography during the last eighty years. This interesting idea has enjoyed a revival through works that use 
extra dimensions to solve the hierarchy problem. A class of these theories use a (4-|-n)-dimensional space- 
time with n compactified dimensions; another class of models considers non-compact extra dimensions. A 
well-known example of this last case is the Randall-Sundrum (RS) model |Q, where the four-dimensional 
world is considered as a three-brane embedded in a five-dimensional bulk. Furthermore in such models a 
four-dimensional graviton exists and is localised within the three-brane. The question which then arises is 
if there is any possibility that other fields, such as gauge fields, fermions and scalars are localised within 
a 3-brane. The problem has been attacked analytically through the search for localized four-dimensional 
fields, where the equations of motion of the bulk fields which are coupled to the background geometry are 
solved p. In this approach, using perturbative tools, a massless photon always appears propagating freely 
in five dimensions. A first attempt of getting evidence of gauge field localization on a brane considering the 
non-perturbative features while using an RS action type for abelian gauge field has been performed in ||3| 
by means of lattice techniques. 

Theories on fiat spacetimes are not, however, devoid of interest. In this paper we continue to explore 
the phase structure of the Abelian Higgs model in five dimensions with anisotropic couplings, defining our 
model on the lattice. In a previous paper ||^ we studied the phase diagram of this model for weak gauge 
coupling in the four-dimensional subspace and found two kinds of layers (or 3-branes) : of Coulomb or Higgs 
type. In the present paper we wish to explore the strongly coupled theory, motivated, in part, by recent 
considerations that indicate that theories studied in ref. |^ are, generically, strongly coupled. 

We find that it is, indeed, possible to tune the lattice couplings in such a way that a phase transition 
from the five-dimensional strong phase to a layered phase occurs. By measuring several order parameters 
we find that the layers are in the Higgs phase and they are separated from each other by a confining force. 
Furthermore we show the precise way in which they are created by studying in detail the phase transition. 
In particular, we find that each layer emerges from the confining phase independently from the others- 
characteristic of a strongly first order phase transition with a very small correlation length in the extra 
direction. 

In section 2 we write down the Abelian Higgs model in five dimensions with anisotropic couplings. In 
section 3 we present the Monte-Carlo results, we exhibit the phase structure and we establish the existence 
of the layered Higgs phase. Also, using mean field techniques, we confirm the Monte-Carlo results and we 
give indication of how the situation changes as we vary the Higgs lattice couplings. 

^ The existence of a layered phase in non-Abehan theories at finite temperature [El has been considered as well. 
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2 Formulation of the model 



The model under study is the Abelian Higgs model in the five-dimensional space. Direction 5 will be singled 
out by couplings that will differ from the corresponding ones in the remaining four directions. 
We proceed with writing down the lattice action of the model. 

S= ^ E (l-cosF^,(x)) + /3;5] (l-cosF^5(a;)) 

X l</j<!/<4 X l</i<4 

+ /?/,5^Re[V(2;)v9(x) - ^ ^* {x)U [,{x)^{x + (i)] 

X 1<M<4 

+ fi'hYl Mv*{xMx) - ^*{x)U^{x)v{x + 5)] 

X 

+ 5][(l-2/3«-4/?,-/?;,V*(x)(^(x)+/3fl(^*(x)(^(x))2], (1) 

X 

where 

F^,^{x) = A^{x) + A^{x + fL)- A^{x + i>) - A^{x), 1 < < f < i, 
Ff,5{x) = A^{x) + A^ix + A) - Af,{x + 5) - ^5(2;) 1 < ^ < 4. 

We have allowed for different couplings in the various directions: the ones pertaining to the fifth direction 
are primed to distinguish them from the "space- like" couplings. The fifth direction will also be called 
"transverse" in the sequel. 

The link variables Up_{x) are defined as e*"s-4s or e^^^rAT respectively, where As, At are the continuum 
fields and as, ax are the lattice spacings in the space-like and the transverse-like dimensions respectively. 
The lattice fields are 

As = asAs, At = otAt- 

In addition, the scalar fields are also written in the polar form (p{x) = p{x)e^^^^\ The order parameters that 
we will use are the following: 

Space — like Plaquette : Ps =< 71:75 E E cosF^,y(2;) > (2) 

X l<fi<u<4 

Transverse — like Plaquette : Pt =< TT75 E E cosF^5(x) > (3) 

X l<ti<4 

Space - like Link : -^^S'=<7^E E cos(x(x + /i) + yl^(x) - > (4) 

X 1<^<4 

Transverse — like Link : Lt ^ cos(x(2; + 5) + ^5(2;) — xi^)) > (5) 

X 

Higgs field measure squared : = ^ P^ix) (6) 

X 

In the above equations is the linear dimension of a symmetric lattice. 

When necessary we will use the order parameters L5 and Ps defined on each space-like volume (layer) 
separately. 

The naive continuum limit of the lattice action |l| may be obtained as follows (where an overbar is used 
for the continuum fields): 

_ haiaT 
^ = 1^^ 



= as\, 1 < < 4, 

Then the transverse-hke field strength 

= Af,{x) + A^ix + A) - Af,{x + 5) - A5{x) (1 < /x < 4) 

goes over to: 

-asiard^Aij^ix)] + aT[asdf,A5{x)] = asar^d^A^ - d^A^). 

Thus 

The space-Uke field strength is treated in a very similar way with the result: 

Flu^a%Fl,, F^, = d^A,-d,A^, l<ii<u<4.. 

This means that the transverse-like part of the pure gauge action is rewritten in the form: 
1 /?o'or 4 r ^2 , 1 (ia'aT 



^ i<^<4 I J ^^^^^ 



On the other hand the space-like part is: 



If we define 



the resulting continuum action reads: 



55 ^T^T 



,2 E -^M^ "'^ „2 E -^m5 



l<//<i/<4 l</i<4 

Defining 7^ = (^)^^^ and using the definitions of f3g, (5g' we find that 

9S as 
I9 = • 

We denote by ^ the important ratio ^ of the two lattice spacings (the correlation anisotropy parameter) 
and finally derive the relation: 

IW 9S . 

After rescaling the scalar fields, one may rewrite the scalar sector of the action in the form: 

2 

/d^^ \D^Tp\^ + ^-^\D,Tp\^+m''r^ + \{Tp*Tp)\ (8) 
1<M<4 ^ 



where D^j^ = 0^ — 1 < /x < 5. 
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We have used the notations: 

7^ 




Ph as Put, 



If we choose a common value for the gauge couphng constants: gs = Qt = 9, (so that 7g = and assume 
that ah the covariant derivatives in equation have the same factor in front: "y^ = the expression does 
not exhibit any anisotropy. However, the naivete of this approach will be manifest by results similar to Ref. 
[P, ^, which indicate that the anisotropy may survive in the continuum limit for a wide range of values of 
lattice parameters both in flat as well as in warped (Randall-Sundrum type) spacetimes. 



3 Monte Carlo Results 

and the Confining— Layered Transition 

For the simulations we use 5-hit Metropolis algotithm for the updating of both the gauge and Higgs fields. 
In order to get better behaviour we use a global radial algorithm and an overrelaxation algorithm for the 
updating of the Higgs field. We simulated the system for 4^, 6^, and 8^ lattices. We made use, mainly, of the 
hysteresis loop method to establish the phase diagram of the system. When necessary, in order to define 
more precisely the phase transition points and study the order of the phase transition we made long runs 
consisting up to 30000 measurements at selected points in the parameter space. 

In whole work we set the four-dimensional gauge coupling fixed at the value Pg = 0.5 and let the gauge 
coupling in the fifth direction run. For some regions of the values of the two gauge couplings we had a 
confining five-dimensional theory. A small value for the Higgs coupling constant in the fifth dimension 
has been chosen = 0.001) (we further discuss this choice in Section 3.2) and we used two values for 
the Higgs self-coupling differing by one order of magnitude: I3r = 0.01 and (3r = 0.1. Thus the phase 
diagram has been found in the f3g — (3^ subspace. 

We study the behaviour of the system in terms of the order parameters defined in Section 2. We proceed 
now with the presentation of the phase structure. 

Figures |^ and ^ show the phase diagrams of the 5-dimensional Abelian Higgs model for the cases of two 
values of coupling, namely Pr = 0.01 and Pr = 0.1, respectively. The phase diagram for both cases has 
been explored in (3'g — (5^ subspace. In Table 1 we show the values for the couplings from which we deduced 
the phase diagrams. The two phase diagrams exhibit similar structure. There are four different phases 
namely the strong phase (S), the Coulomb phase in five dimensions (C5) and the Higgs phases in four (^^4) 
and five dimensions {H^). More details about the nature of the different phase transitions will be given 
later. The crucial feature, however, is the existence of a phase transition between the strong phase (eg. a 
confined phase in five dimensions) and the i?4 phase -a phase of broken U(l) symmetry in four dimensions- 
which gives rise to a constitution of a layered phase with broken symmetry on each layer. Let us give now 
some representative results which lead to the identification of the various phases of the model. 

• The S — H4 and S — H5 phase transitions 

In Figure ^(a) , we present the hysteresis loops concerning the space-like plaquettes, as f3h runs, for (3r = 0.01 
and for two values j3'g = 0.2 and fi'g = 0.7. In Figure ^ (b) the corresponding results for the transverse-like 
plaquettes are depicted. 
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Figure 1: Phase diagram for 
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Figure 2: Phase diagram for = 0.1 



The behaviour of the space-hke plaquettes and also the behaviour of R^, shown in Figure lead to the 
conclusion of a phase transition between the five dimensional strong phase and a phase with broken symmetry 
which probably is of first order by virtue of the large hysteresis loop. In addition, the transverse-like 
plaquette, Pt, for P'g = 0.2 remains almost constant to a small value (it equals the value Pg/2, labeling the 
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Pr = 0.01 


= 0.1 






S — i?4 


0.05, 0.10, 0.20, 
0.40, 0.45 


0.10, 0.20, 0.25, 
0.30, 0.40, 0.45, 0.46 


S-H5 


0.50, 0.70, 0.80, 
0.90, 1.00 


0.47, 0.50, 0.70, 
0.90, 1.00, 1.10 


C5 — 


1.40, 1.50, 1.90 


1.20, 1.30, 1.50 


H/^ — 


0.40, 0.50, 0.60 


0.50, 0.60, 0.80 


S-C5 


0.05, 0.15, 0.20, 
0.25, 0.30 


0.10, 0.20, 0.30 



Table 1: The phase which the phase diagrams given in Figures land I are based on. 
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0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 

P. P. 
(a) (b) 



Figure 3: Hysteresis loop for space-like (a) and transverse-like (b) plaquette for two values of (3g = 0.2 and 
0.7. 



Strong phase) while the corresponding one for f3g = 0.7 increases with f3fi as a phase transition occurs. This 
is a serious indication that there are two different Higgs phases: in particular one is a Higgs phase in four 
dimensions (with confining behaviour along the fifth dimension) and the other is a five-dimensional Higgs 
phase. 

• The H4 — H5 phase transition 

The fact that the fifth dimension is confining is made clearer from the result of Figure ^: we keep f3h = 0.40 
and let P'g run. The results depicted in this figure correspond to the values of transverse-like link and 
transverse-like plaquette which are small enough for small values of P'g and they increase as this coupling 
parameter is running to larger values. It can be noticed that there is no obvious hysteresis loop formed and 
the system passes from the H4 (where the two order parameters take small values) to the phase in a 
fairly smooth way. 
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Figure 4: Hysteresis loop for R^, j3'g = 0.2,0.7. 



V=6' Pk=o.oi 
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0.0 -I 1 , 1 , 1 , 1 , 1 , 

0.1 0.2 0.3 0.4 0.5 0.6 



K 

Figure 5: Hysteresis loop analysis for Lt and Pt in the region between and phases. The transition 
is not showing formation of a loop and seems to be fairly smooth. 

• The S — C5 and C5 — H5 phase transitions 

The existence of the C5 phase is indicated in Figure ^ which contains the hysteresis loop for Ps and Pt for 
running (3'g. As it can be easily seen, these two quantities pass from a region where their values are almost 
half of their corresponding gauge couplings (Strong phase) to a phase (C5)where their values tend to one. 
The large hysteresis loops indicate a first order phase transition for S — C^. 

Finally, the transition between C5 and is shown in Figure |^ where we set (3'g to 1.5 and let (3h run. 
In Figure |^ (a) we show the behaviour of the space-like and the transverse-like link from which we can 
see that they both exhibit a gradual increase as j3h grows. In addition, the corresponding behaviour of iil^ 
values in Figure ^ (b) implies the transition from a five-dimensional Coulomb phase to a five-dimensional 
Higgs phase. 

The order of the phase transition is not obvious and further study will be needed. To this end we measure 
the susceptibilities for the Ls order parameter at the value fi'g = 1.5. The results shown in Figure ^ seem to 
be consistent with a second order phase transition as the maximum values of susceptibility for each volume 
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(a) (b) 
Figure 7: Hysteresis loop for space-like and transverse-like link (a) and (b), for /?' = 1.5. 



grows with some power of the volume which is obviously much less than unity. 

Up to now, we have given some examples of the behaviour of the system for = 0.01. In general, 
the phase structure for I3r an order of magnitude larger than this , e.g. Pr = 0.1, is similar. However, we 
can elaborate on some points concerning the structure of the layered phase which is formed in the S — H4, 
transition on one side and the S — transition on the other. 

In Figure ^ we can see an example for the transition from S to i/4 {fi'g = 0.2) and from S to {j3'g = 0.9), 
by giving the hysteresis loops for Ps and Pt- The difference of this figure with Figure ^ consists in the 
"weaker" transition from S to Hj^ phase, since the loop is much smaller. A more detailed study for the case 
I3'g = 0.2 shows (Figure [l0| ) that a clear hysteresis loop is formed, indicating a first order phase transition, 
though seemingly weaker, than that of 13r = 0.01 case. However, the trasverse-like order parameter Pt does 
not change at all as (3^ increases. This fact points out that a layered phase has appeared and the layers are 
decoupled from each other. 
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0.280 0.282 0.284 0.286 0.288 0.290 0.292 



Figure 8: Susceptibility S{Ls) for three values of lattice volume for P'g = 1.5 (the error-bars are smaller 
than the sizes of the symbols). 



V=6^ P^=0.1 



V=6= P^=0.1 





0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 



(b) 

Figure 9: Hysteresis loop for space-like (a) and transverse-like plaquette (b) for f3g = 0.2 and f3g = 0.9. 
3.1 Multi-layer Structure 

The next step is to study the behaviour of the layers one by one as the system moves through the phase 



transition. This is shown in Figures 11 (a) and (b) which correspond to S — H/^^ and S — phase transitions 
respectively. In the Figure ^ (a), corresponding to the S — case, one may easily see a "non-sychronised" 
transition exhibited by L5 defined on each space-like volume (layer) in contrast with Figure ^ (b), in which 
the corresponding order parameters indicate the phase transition simultaneously (Obviously, in figure (b) 
the hysteresis loops formed by the Lg defined on each layer can not be distinguished from the corresponding 
hysteresis loop due to L5 defined on the volume.) This specific behaviour of hysteresis loops, may actually 
serve as a "criterion" to characterise the layered phase. 

We also reach the same conclusion by considering the [5r = 0.01, [3'g = 0.2 case, for a 4^ lattice in more 
detail. The result, shown in Figure leads to the same conclusion: The very existence of the layered 
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p„-0.1 

0.9-1 , 1 , 1 , 1 , 1 — ^ r- 




0.2 -| , 1 , 1 , 1 , 1 ^ 1 , 1 , 

0.440 0.445 0.450 0.455 0.460 0.465 0.470 0.475 

Figure 10: Hysteresis loop for space-like link and plaquette for = 0.2 in more detail, indicating a first 
order phase transition. 



V=6" P^=0.1 6=0.2 




V=6 13^=0.1 13-0.9 




0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50 



(a) (b) 

Figure 11: Hysteresis loop for P'g = 0.2 (a) and f3g = 0.9 (b). We show the very different way for the 
transition of Ls calculated either on the five dimensional volume and on four layers. In the figure on the left 
the S — H4 transition is presented where the layers show a decoherent behaviour on the phase transition in 
contrast with the S — phase transition on the right where the transition for Ls is identical for all layers, 
so it coincides with the mean values over the 5D-volume. 

phase (and consequently of confinement in the fifth dimension) shows up in the "incoherent" behaviour of 
the space-like volumes (layers) as the phase transition takes place. 

This behaviour can be further confirmed by long runs in the transition region. It is to be expected that 
the S — Hi phase transition is of first order. Therefore, we would expect a two peak signal in the order 



parameter disrtibution at equilibrium. Nevertheless the situation, shown in Figure 13 (a) concerning the 
distribution for the order parameter Ls for a value of f5h near the transition region for 1/ = 6^ is by no means 
what one would expect normally. The multipeak structure seems rather strange. It should be noticed that 
the same occurs for the order parameter Ps too. However, the study of the same order parameter defined 
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1.0 - 




0.0 -I 1 1 1 1 1 1 1 1 1 1 1 1 

0.28 0.30 0.32 0.34 0.36 0.38 0.40 

P, 

Figure 12: Hysteresis loop for Ls showing the way that the layers decouple for (Sr = 0.01 at P'g = 0.2. 

on each space-like volume is more illuminating. For example in Figure |l^ (b), we can see the distribution 
of Ls values on each layer. We show four out of six distributions of Ls corresponding to the four space-like 
layers within the five-dimensional volume. The distributions which are produced appear quite usual and 
they show that at the same time one layer is in the strong phase (called 2nd in tlie figure), another has 
already passed to the broken phase (3^^) and others produce a two peak signal result. The result is that 
the strange picture of the distribution formed in Figure ^ (a) is resolved if we analyze the behaviour of the 
system on each layer as the system undergoes the phase transition. 




(a) (b) 

Figure 13: Distribution for Ls on the volume (a) and on each layer (only four of the six are shown in the 
figure) (b) for 6^ lattice in the critical region for the S — phase transition. 

This result is found for all of the volume sizes (e.g. 4^,6^,8^) which we have worked on. Although it 
is consistent with a first order phase transition it lends support to the view of a decoherent behaviour for 
every four-dimensional volume (layer) in the transition region between five-dimensonal strong phase and 
the four-dimensional layered Higgs phase. This certain behaviour describes a dynamical decoupling of the 
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layers and provides a possible mechanism for localization of the fields on the layers. 



3.2 Mean Field Approach 

The Mean Field Approach provides a point of view complementary to Monte Carlo sumulations. Although it 
is, by construction, blind to spatial fluctuations, and cannot, therefore, see the multi-layer structure, it does 
not suffer from the finite volume effects, that limit Monte Carlo simulations. Furthermore, it is expected to 
be a reliable guide to the phase diagram, the higher the dimension of the system under study. 

Thus, in this work we use the Mean Field analysis (i) to show that the small value of we have chosen is 
really a suitable one to reveal the layered Higgs phase and (ii) to provide evidence that as Higgs self-coupling 
takes larger values -up to = 0.2- the phase transition from the Strong to the H4 phase may become 
second order. 

We start with the action (jl]). We fix the gauge by imposing U^{x) = I and use the translation- invariant 
Ansatz ||, ^, Ujx{x) = w, 1 < ^ < 3; U'^{x) = v' . We also introduce the variables for the Higgs field, 

(t>{x) = p(x)e'>^(^) (9) 

and have also assumed a translationally invariant Ansatz^ p = p{x), x = xi^)- The free energy, which should 
be minimized to get the mean field solution, reads: 



2. ,2 
X 



+(1 - 2/?fi)p2 + _ hiog[p'] 
+ "iava - 3/05 [Jo (a)] + a'v'^, - log[lQ{a')] + x^x " ^og[Io{x)] (10) 

The parameters a, a' and x are conjugate to Va, v'^, and respectively. 

We study the phase structure using the order parameters defined in Section 2. 

Our Monte Carlo analysis has been performed setting = 0.001. It is of interest to have an idea what 
happens when takes other values. 



In Figure |14| we depict the behaviour of Pt setting /3/j = 0.01 and for the two cases corresponding to the 
transitions from the Strong phase to and H^. The corresponding values for are 0.6 and 0.2 respectively 
(see, also, Figure |l]). From Figure ^ (a) we can see that as increases the transition behaviour to 
is fairly the same. On the contrary. Figure ^ (b) shows that the increase in leads to a big increase 
in the value of Pt which becomes compatible with the value characterizing the H^, phase. This provides 
an indication that at some value, (3'^ ~ 0.1, the phase transforms to H^: as (3'^ increases extends, 
covering the region occupied before by 



Figure |15| provides evidence of the expected weakening of the Strong - phase transition as (5r 
increases. It reveals that, although for I3ii = 0.01, Ps exhibits the same behaviour for and phases, 
for f3ii = 0.2 the picture changes substantially. The S — phase transition provides a sign of a smoother 
phase transition leading probably to second order. Mone Carlo simulations in this region of parameter space 
are required to complement this evidence. 
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Figure 14: We show the behaviour of Pt as the value of grows. The transition from Strong phase to the 
i?5 phase does not seem to change (a), in contrast with the one to from Strong to which is washed away 
(there is a radical change to the value of Pt at a critical point of f3h) leading thus to the existence of only 
a 5D-Higgs phase. 
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P, = 0.2 p'=0.6 



Figure 15: As the value of increases the S — phase transition becomes weaker, as expected, while the 
S — transition remains first order. 



4 Conclusions 

In this paper, we have shown, using Monte-Carlo methods, that a layered Higgs phase actually exists in the 
phase diagram of the strongly coupled five-dimensional Abelian Higgs model and that it emerges from the 
confining bulk phase through a first order phase transition. In fact we find the existence of multi-layers, as 
each brane passes from the confining phase to the Higgs phase independently of the others. 

Using mean field techniques we performed a scan of the phase diagram as the Higgs couplings changed 
and found evidence that, as the Higgs self-coupling increases, the emergence of the Higgs layers from the 
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confining bulk phase softens and may become second order, leading to new continuum theories. Further 
work will clarify this issue. 

It is worth stressing that the layers of our model are, indeed, 3-branes. In string theory one expects 
symmetry enhancement when branes coincide and the question arises, whether such an effect could be visible 
within our field theory context. This effect is the manifestation of new, non-perturbative, degrees of freedom. 
In our context this would mean introducing magnetic monopoles, that would promote these 3-branes into 
true D-3-branes. One way of achieving this could be using twisted boundary conditions, along the lines of 
ref. i. 

This raises the question of what aspects of our study may also be of relevance to Yang-Mills theories ||^, 
10|. As is well known, these are confining in less than five dimensions, so do not seem to admit four 
dimensional layered phases. The simplest example would, thus, be an anisotropic SU{2) theory in six 
dimensions and in this case the layers form an unnatural five-dimensional Coulomb phase . The interesting 
point would be if the SU(2) Higgs model with extra dimensions could reveal a layered Higgs phase in four 
dimensions which bears a closer resemblance to the SM. Monte Carlo simulations of such theories directly 
are, unfortunately, inconclusive with current technology-though the use of massively parallel clusters may 
offer some hope for the not too distant future. Using suitably reduced models, on the other hand, such as 
the one studied here, could provide some useful information on their structure. 
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